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APPLICATION OF NEMAN TRANSFER FUNCTIONS TO A DESIGN PROBLEM 


By James J. Adams 
NASA Langley Research Center 

SUMiyiARY 

An analytical design study was made of a proposed full-scale^ manually 
controlled lunar landing simulator using analytical transfer f\mctions for the 
pilot control response along with the analytical representation for the mech- 
anisms. The simulator reproduced the lunar environment by supporting five- 
sixths of the weight of the test vehicle with an overhead cable. The cable was 
kept directly over the test vehicle by the automatic control of the longitudinal 
drive mechanism of the simulator. The results showed that the dynamic charac- 
teristics of the simulator that could be expected in the actual system were in 
a range that would influence the response of the manually controlled systems 
which were to be tested. 

When the simulator was put in operation^ the results of the analytical 
study were checked. The simulator was operated with the gain of the longitudi- 
nal drive set as high as was feasible with the actual mechanism and with a low 
gain to determine if this change would affect the pilot's response. The pilots 
reported that the degraded system was more difficult to control^ and the records 
clearly showed a decrease in system damping with the degraded system. 

INTRODUCTION 

One of the reasons for determining human transfer functions is so that 
evaluation and prediction of the performance of manually controlled systems 
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can Tdc accomplished during design studies. This report explains how such 
transfer functions were used in the design analysis of the drive system of a 
lunar landing simulator. 

The simulator was designed to provide a i+OO-foot hy 50-foot hy l80-foot 
high volume in which lunar landing maneuvers co-uld he studied. G?he lunar 
gravity was simulated hy supporting five-sixths of the weight of the test 
vehicle hy a cable. The load in the cable was regulated hy measuring the load 
with a strain gage and operating the overhead winch in response to the error 
in this measured load. The cable was kept directly over the vehicle hy meas- 
uring the cable angle at the overhead, traveling bridge, and moving the bridge 
in response to this measured angle. It is this longitudinal drive system of 
the bridge that is the subject of the design study reported in this paper. 

The transfer functions used to describe the pilot's control action were 
derived in reference 1. These transfer functions describe the pilot's control 
used when controlling the multiloop system representative of the lunar landing 
horizontal translation maneuver. The use of these pilot transfer functions in 
determining the most suitable simulator drive characteristics will be presented. 

SYMBOLS 

mass, slugs 
translation, ft (m) 
pendulum length, ft (m) 
pendulum angle, deg 

gravity, 52.2 ft/sec^ (9«8l m/sec2) 
simulator drive system control gains 



s 


LaPlace operator, per second 



e 


pitch attitude angle, deg 
S control moment, rad/ sec^ 

u) undamped natural frequency, rad/sec 

h maximum deflection of first cable vibration mode, ft (m) 

Ki,K 2 ,t gains in analytical transfer function of pilot 
leakage flow, in5/sec 
relief valve flow, in5/sec 
V oil volTime under compression, in5 

B bulk modulus of oil, psi 

Pjjj^ motor pressure, psi 

Djji motor volume, in5 

Jjn motor Inertia, in-lb -sec^ 

motor rotation, rad/sec 
n gear ratio 

Bjjj motor damping, in-lb/rad/sec 

Kl tire torsional spring constant, in-lb/rad 

Tl tire torque, lb 

Jjj load Inertia, in-lb -sec^ 

Bjj load visions friction, in-lb/rad/sec 

Subscripts : 

B bridge 

y vehicle 

c command 

e error 
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DESCIPTIOK OF PROBLEM 


Simulator 

Simplified equations of motion for the longitudinal motion of the over- 
head bridge, and pendulum consisting of the cable -supported veh i cle are: 

(me + “v)% + ^ 

myl% + 10^1^9 + Ittyglcp = 0 

!Ehe characteristic equation for the system in LaPlace rotation is 

(wq + i^v) + (“b + ™v)(“7Sl) - ° 

which' is the familiar equation for a pendulum in which the frequency is pri- 
marily determined by the length 1 but with an additional term which repre- 
sents the influence of the bridge being free to move, which reduces the fre- 
quency somewhat. 

If the control commands that the bridge accelerates as a function of cp 
and 9 then the equations of motion become: 

(®B + 

mylitg + myl cp + IQyglCp = 0 
The characteristic equation, in LaPlace rotation, is: 

(mg + my)^myl2js2 + + ^mg + my^^mygl) - my2l2s2 + m^lK^ = 0 

It can be seen from this equation that the gain, in the coefficient of 

of drive system control will supply damping to the system, and the Kq> gain 



will increase the penduliom frequency so as to keep the bridge above the sus- 
pended vehicle. 

The bridge drive unit was an electro-hydraulic unit consisting of a syn- 
chronous electric motor which drove a variable displacement hydraulic pump. 
Control of the bridge was exercised by the operation of a pump stroker, which 
controlled the displacement of this pump. The p\imp drove the fixed displace- 
ment hydraulic motors attached to the wheels of the bridge. A fixed displace- 
ment of the stroker produced a steady state constant velocity of the bridge. 

Since the stroker controlled bridge velocity instead of bridge accelera- 
tion, the control function of accelerating the bridge as a function of pendulum 
angle was achieved by displacing the stroker as a function of the Integral of 
the pendulum angle, and the fuinction of accelerating the bridge as a function 
of rate of change of pendulum angle was achieved by displacing the stroker as 
a function of pendulum angle. 

The predominant dynamic characteristic of the drive unit was the oscilla- 
tory response of bridge velocity to stroker displacement that resulted from the 
compressibility of the hydraulic fluid as it reacted against the mass of the 
bridge. This dynamic characteristic is expressed by the equation 

__ ^ 


Stroker displacement g 




,2 + + 1 


where 

=1 

^m 

B 

V 


oil volume delivered by pump per unit displacement of stroker, in3/unit 

motor displacement, 28 in3/rad 

bulk modulus of oil, 1 X 105 psi 

oil volume under compression, 600 in3 
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Jjj total inertia of load reflected at motor output shafts 5050 in-lb/sec^ 

Kl leakage coefficient, in5/sec/psi 

For the system under study, the natural frequency of the drive unit as deter- 
mined hy this equation is 4.6 rad/sec. The detailed analysis of the system 
naturally included this characteristic. 

In addition to the drive unit dynamics, the following dynamic factors were 
also Included in the detailed analysis 

(1) the compliance of the primary electric motor 

( 2 ) the time constant and the limit displacement of the stroker 

( 5 ) the motor leakage, as a function of ^pressure 

(4) pressure relief valve flow 

( 5 ) nonlinear friction of the gearing 

(6) tire compliance 

Another important dynamic characteristic of the system was the oscillatory 
characteristic of the cable. Cable vibrations added to the measured cable angle 
and therefore added a spurious signal to the control signal. Equations for the 
first two modes of vibration for different, fixed cable length were determined 
and used in the analysis. These equations Included the effect of a lumped mass 
located near the vehicle which represented the whlffletree which was a part of 
the support and gimbal arrangement of the vehicle. The frequency of the first 
mode for a 200-foot cable length was 8.72 rad/sec, which is very close to the 
natural frequency of the drive unit, and which therefore put a limit on the pre- 
cision with which the bridge could be maintained above the vehicle. This vibra- 
tion frequency would Increase at shorter cable lengths, and also change with 
vehicle wei^t. The analysis was made for fixed cable lengths and vehicle 
weights. 
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Pilot Transfer Function 


The pilot transfer functions^ which were used in conjunction with the 
simulator equations, are derived in reference 1 and are repeated here. The 
pilot -vehicle system involved in the landing maneuver is a multiloop system 
described in the block diagram presented in figure 1. The inner loop deals 
with the attitude control of the vehicle. The vehicle response to attitude 
control was assumed to contain a proportional rate feedback, and is given by 
the equation 

e ^ o.b 

S s(s + 0 . 5 ) 

which defines a rate system with a rate response with ,a time constant of 
2 seconds. Eeference 1 demonstrates that a pilot's response in such an inner 
loop is given by 

5 ^ ^ 96(1 + 0.4s) 

®c,e (s + x)^ (s + 6)^ 

The combination of the pilot and vehicle gives this inner loop a closed-loop 
characteristic frequency of 1.2 rad/sec and a damping ratio of 0 . 26 . 

The outer loop of the system deals with the longitudinal translation, and 
the vehicle response to attitude angle is given by a pure inertial response 

X _ 36 

e " s2 

This relation is derived from the linearized equation of motion for the hori- 
zontal component of acceleration due to one-sixth of thrust that would be in 
effect in the lunar environment 
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X = ^ sin 0 
6 


using small angle linearization 


X = ^9 = 5*369 ft/sec2 = 1.630 m/sec2 
6 


X 


^=5.36 


ft/ 


sec‘^ 


rad 


= 1.63 


m/ sec^ 
rad 


Reference 1 demonstrates that the pilot response in such an outer loop is 

0c ^ 0.9(1 + 9*2s) 

(s + 10)2 

This pilot response defines a characteristic response of the complete system 
■which has two small real roots, s = -O.167 and s = -0.336* In terms of an 
oscillatory response/ these two roots define an overdamped response with a 
natural frequency given hy 

m = J(0. 167) (0.336) = 0.236 rad/sec 

This system frequency characterizes the translation response of the system. 

Since the response characteristic of the longitudinal drive system of the 
simulator must have a response frequency hi^er than that of the system which 
is to he tested, this calculated response characteristic of the pilot controlled 
translation response provides a first, rough criterion for the required charac- 
teristics of the longitudinal drive system. 

ANALYSIS AND TEST RESULTS 

A detailed analytical study was made to determine precisely the drive 
system characteristics and the suitahility of the characteristics. The computer 
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diagram used in this study is presented in figure 2 . It includes the complete 
representation of the drive system^ the cable dynamics, and the pilot controlled 
vehicle representation. 

The first phase of the study was conducted to determine just how fast the 
bridge could be made to respond. Open -loop step thrust inputs to the vehicle 
were used as the forcing function in these studies. The results showed that 
the presence of the cable vibration modes of motion in the system placed an 
upper limit on the pendulum damping gain K^. If this gain was adjusted too 
high, the first vibration mode would become unstable, as is illustrated in 
figure 3. 

The limit on the pendulum danp)lng gain placed further restriction on the 
pendul-um frequency gain Kqj. The range of possible system characteristics that 
could be achieved is shown in figure h, which shows the vehicle velocity response 
to a 2 -second thrust impulse. The oscillatory natiire of these responses is the 
result of the bridge drive system characteristics. It can be seen that a well- 
damped response with a frequency of 1-57 radlans/sec (a period of 4 seconds), or 
a poorly damped response with a frequency of 2.5 rad/sec (a period of 2.5 sec- 
onds) could be achieved. 

Both of these frequencies are above the O.236 rad/ sec frequency for the 
pilot controlled translation response of the lunar landing system. However, it 
cannot be confidently concluded that they are sufficiently hi^ so as to have 
no effect on the simulation. To determine what effect the bridge response 
might have on the pilot -controlled maneuver, the analytical representation of 
the pilot and vehicle were included in a closed-loop representation of the 
complete system, and the response to a commanded 200-foot displacement was 
deteimlned. The results are presented in figure 5 ^ which show, first,- the 
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response of the pilot-vehicle combination alone ^ and then the response of the 
pilot-vehicle -simulator combination with the two different simulator character- 
istics presented before. Using the pilot-vehicle conbination response as the 
standard for comparison it can be seen that including the simulator bridge 
dynamics in' the loop does Indeed influence the response. With the lower gain 
bridge control the system is degraded to the point of instability. It was 
therefore concluded that the simulator should be adjusted so as to have as high 
a frequency characteristic as possible. Also, it was indicated by the analysis 
that the pilots might find the simulator slightly more difficult to control 
than the real lunar landing system. 

Post-Analysis Tests 

When the simulator was put in operation, the system characteristics which 
could be achieved with the actual mechanism were determined. It was found that 
the highest stable pendxilm frequency that could be obtained was 1.4 rad/sec 
(a period of 4.5 seconds). Piloted i*uns were then made with this highest fre- 
quency response of the bridge, and with the Kqj gain placed at a lower setting 
to check the analytical results that such a change would affect the piloted 
system response. The lower frequency used was approximately 0.8 rad/sec 
(a period of 8 seconds) . Flg'ure 6 shows cable angle responses of the simulator 
in these two conditions to open-loop step inq)ulses. The piloted tests were 
started with the vehicle hovering at an altitude of approximately 50 feet. The 
pilot then translated the vehicle 200 feet and attempted to stop and hover over 
a mark located on the ground. The run with the higher response was made first, 
and then immediately a second run was made with the lower setting. Two differ- 
ent pilots were used. 
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Sample time histories of these tests are shown in figure 7- In this test 
the first maneuver, from the 50 -foot point to the 275 -foot point, was done with 
the high-gain drive system. At the 50-second mark, while the pilot was turning 
180 °, the drive system gain was readjusted to the low-gain setting, and the 
pilot started to go hack to the 50-foot mark. It can he seen that the general 
nature of these maneuvers is very similar to that computed in the analytical 
study. With the lower response characteristics for the simulator a very 
noticeable decrease in damping of the attitude angle can he seen. The pilot 
did not continue the maneuver in this case, hut rather dropped the intention 
to precisely control translation. He stopped the attitude oscillation and 
landed at the point that was helow him at that time. The pilots commented that 
in the lower response riins they were having more difficulty in controlling the 
vehicle, and as one pilot said, he felt he was in a "pilot-induced oscillation" 
condition. 

Since the simulator characteristics that were achieved with the actual 
mechanism when it was put in operation were not the same as those determined in 
the analytical study, and since the pilot's response showed a lower attitude 
angle limit in the flight tests than was assumed in the analytical study, the 
analytical study was repeated in an attempt to reproduce more closely the fll^t 
time histories. In these repeated calculations the drive system gains were 
adjusted to give a poorly damped 4.5-second pendulum period in one case and a 
well damped 8 -second pendulum period in the second case to correspond to the two 
conditions that were tested in the flight tests. The same linear transfer func- 
tions were used for the representation of the pilot. The attitude limit was 
placed at 10 °, which corresponds more closely with the limit used hy the pilot in 
the flight tests than did the 40° limit used in the Initial analytical study. 
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With this 10° limit on attitude angle^ a ■well- controlled attitude time history 
was calculated^ shown in fig-ure 8 (a) ^ which is a better reproduction of the 
flight time history than was the initially calculated response. When the simu- 
lator characteristics were changed so as to have an 8 -second period, a deterio- 
ration in the stability of the calculated attitude angle time history resulted. 
These repeated calculations further confirm the conclusion drawn from the 
initial analytical study. 

CONCLUSIONS 

The experience gained in the exercise of predicting the characteristics of 
a manually controlled simulator and verifying these characteristics when the 
simulator was put in operation has demonstrated the -validity and usefulness of 
analytical expression of human response. The analysis showed that the range of 
simulator dynamics which was likely to occur in the simulator would influence 
the response of the pilot -controlled maneuver, and tests with the completed 
hardware confirmed this conclusion. 

It was concluded that the gain of the simulator longitudinal drive system 
should be kept as high as possible to minimize the effect on the piloted 
maneuvers, and it was indicated by the analysis that the tasks performed with 
the simulator might be sli^tly more difficiolt than the same tasks performed 
in the 1 -unar environment. 
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Figure 4.- Calculated response of simulator to a 2-second thrust impulse. 



(a) Pilot vehicle combination. 



(b) Pilot-vehicle-simulator combination with 
high gain simulator. 



(c) Pilot-vehicle-simulator combination with 
low gain simulator . 


Figure 5.- Calculated vehicle attitude response to a 200-foot step 

translation command. 
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Pilot-vehicle- simulator combination with low gain 
(4 second period) simulator. 



Time, sec 

(b) Pilot- vehicle- simulator combination with very low 
gain (8 second period) simulator. 


Figure 8.- Repeated calculation of vehicle attitude response to a 
200-foot step translation command. 



